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CO , Abstract. We present a comprehensive study of the z = 0.43 host galaxy of GRB 990712, involving ground-based photome- 

try, spectroscopy, and HST imaging. The broad-band UBVRIJHKs photometry is used to determine the global spectral energy 
distribution (SED) of the host galaxy. Comparison with that of known galaxy types shows that the host is similar to a moderately 
reddened starburst galaxy with a young stellar population. The estimated internal extinction in the host is Ay ~ 0.15 ± 0.1 
and the star-formation rate (SFR) from the UV continuum is 1.3±0.3 Mq yr^^ (not corrected for the effects of extinction). 
, Other galaxy template spectra than starbursts failed to reproduce the observed SED. We also present VUT spectra leading to the 

detection of Ha from the GRB host galaxy. A SFR of 2.8±0.7 Mq yr~^ is inferred from the Ha line flux, and the presence of 
a young stellar population is supported by a large equivalent width. Images from HST/STIS show that the host has two separate 
knots, which could be two distinct star-forming regions. 
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K*" . 1. Introduction to noise ratio spectra. Ground-based broad-band photometry 

' presents a useful alternative for investigating the spectral en- 

For all but one of the Gamma-Ray Bursts (GRBs) where the ^^.g^ (jistxibution (SED) 



o 

cn 



^ position of the X-ray, optical, or radio afterglow has been lo- 



calised to an accuracy of less than 1", follow-up deep obser- 
vations have revealed underlying galaxies. The current sam- The relatively small size of GRB hosts sometimes makes 
pie consists of --45 such GRB hosts, and --35 of these it difficult to tell which morphological type it is, e.g. whether 
have measured redshifts in the range 0.168 < z < 4.5 the radial intensity profile of the galaxy is best fitted by an ex- 
(Hjorth et al. 2003; Andersen et al. 2000) and have magnitudes ponential disk profile or an elliptical profile. However, studies 
21 < R < 30. The faintness of the hosts requires long in- of some hosts with the HST have shown that an exponential 
tegration times on the largest telescopes to obtain high signal profile provides a good fit to the surface intensity distribution 

(Odewahn et al. 1998; Fruchter et al. 2000; Hjorth et al. 2002). 

* Based on observations collected at the European Southern Comparing the SED of the host with the SEDs of known galaxy 

Observatory, La SiUa, with the ESO 3.6m telescope, NTT and VLT ^ provides an alternate method of estimating the galaxy 

(ESO Programmes 165.H-0464(I) and 165.H-0464(E), and 265.D- ^^p^ g^j^^j^^ ^2001) analysed 6 GRB hosts in this way, 

5742(B)) by the Gamma-Ray Burst Afterglow Collaboration at ESO , . , n ^ i i i r,r-i^ i • ^ i 

^ , , , , ■ ...x^ .,,^ showing that all of them had SEDs characteristic of starburst 

(GRACE) group. Further based on observations with the Danish 1.5m . ° , . , ,t^t^^ ■ , ■, 

^ , galaxies. This is expected it GRBs are associated with massive 

tClCSCODC 

Based on observations with the NASA/ESA Hubble Space collapsing stars as suggested by Woosley (1993) and Wijers 
Telescope, obtained at the Space Telescope Science Institute, which etal. (1998) andrecendy observed for the GRB 030329 (Stanek 
is operated by the Association of Universities for Research in et al. 2003; Hjorth et al. 2003). A small age for the burst pop- 
Astronomy, Inc. under NASA contract NAS5-26555. ulation gives an indication that GRB progenitors ai'e massive 
Correspondence to: lchristensen@aip.de Stars, whereas SED ages much longer than the fife times of the 
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most massive stars could indicate a binary merging event as the 
cause of the GRB (Eichler et al. 1989). 

The GRB 990712 host is bright relative to other GRB hosts, 
and therefore serves as a good case for studying the multi- 
band SED. This paper is one in a series of papers on the 
SEDs of GRB hosts. Studies of the hosts of GRB 000210 and 
GRB 000418 are presented in Gorosabel et al. (2003a) and 
Gorosabel et al. (2003b), respectively. 

The previous studies of the GRB 990712 afterglow and the 
host are summarised in Sect. 2. In Sect. 3 and Sect. 4 we present 
photometry and spectroscopy of the GRB 990712 host. The 
morphology of the host is investigated in Sect. 5. A compari- 
son of the SED derived from all the observations with spectral 
synthetic templates is described in Sect. 6. In Sect. 7 we es- 
timate the SFR of the host galaxy using two SFR estimators; 
first using the UV continuum and second the Ha hne flux. In 
Sect. 8 we discuss the results. 

Throughout the paper we assume ilm = 0.3, JIa = 0.7 and 
Ho = 65kms-iMpc-\Atflieredshiftoffliehost,0 = 0.433, 
the luminosity distance is di, = 7.93 x 10^''cm. 

2. GRB 990712 

GRB 990712 was detected on July 12 1999 at UT 16:43:02, 
by the Italian-Dutch satelHte BeppoSAX (Frontera 1999). The 
burst had the strongest afterglow observed in X-rays to that 
date. Its afterglow was found approximately 4 hours after the 
trigger by Bakos et al. (1999). A spectrum obtained shortly af- 
ter revealed a redshift of z = 0.433 from the emission lines 
[O II], [O III], H7 and H/3 as well as Mg I and Mg II absorption 
lines (Galama et al. 1999; Hjorfli et al. 2000). The relatively 
low redshift makes it one of the closest GRB hosts. Only the 
GRB 980425 at ^ = 0.0085 (Galama etal. 1998), GRB 011121 
atz = 0.36 (Infante et al. 2001) and GRB 030329 atz = 0.168 
(Greiner et al. 2003) were nearer. 

A spectrum of the combined flux from the host and the af- 
terglow was obtained by Vreeswijk et al. (2001a) 1 .5 days after 
the burst. Because of the brightness of the host, the spectrum 
shows distinct absorption lines and emission lines from the host 
itself. 

The [Oil] emission line flux was measured to be 
(3.37±0.2) X 10-16 erg cm'^ s'^ (Vreeswijk et al. 2001a), 
which corresponds to a SFR of 2.7±0.8 MQyr^^ using the con- 
version from measured flux to a SFR from Kennicutt (1998) 
(Hereafter K98). Converting the flux at restframe 2800 A 
to a SFR gives a similar result. An internal extinction of 
was inferred from the flux ratio of the hy- 
drogen lines H7/H/3. The host was observed at radio frequen- 
cies (1.4 GHz) by Vreeswijk et al. (2001b), who did not find 
any radiation from the host to a limit of 70 /xJy. This upper 
limit implies that the total unextincted SFR in the host is less 
than 100 Mq yr"^. This is in great contrast to the measure- 
ments of the host of GRB 980703, which was found to have 
a SFR of ~500 MQyr"^ measured from its radio flux (Berger 
et al. 2001; Berger et al. 2003). The GRB 990712 host has an 
IR-luminosity which is 20 times less than that of the luminous 
GRB 980703 host. This suggests that different types of galaxies 



can host GRBs some having more dust enshrouded star forma- 
tion than others. 

3. Imaging and photometry 

We have analysed both ground based and HST images of the 
host of GRB 990712. The ground-based observations consist 
of UBVRIJHKs images obtained at different dates and using 
different instruments. All of the data presented here were ob- 
tained more than one year after the burst, so that the flux con- 
tribution from the afterglow is negligible. The data obtained 
from the Danish 1.5-m in September 2000 consist of images 
in Bessel B, V and R and Gunn I filters, and the U band data 
were obtained at the ESO 3.6-m telescope the night of Aug. 
13, 2001 using the EFOSC2 instrument. Near-IR JHKs images 
were obtained at the NTT with the SOFI instrument over two 
nights from Aug. 1-2, 2001 . Only the first night was photomet- 
ric according to the ESO webpage'. 

In order to obtain a more reliable optical estimate of the 
SFR of the host of GRB 990712, we performed spectroscopic 
observations centered on Ha, using FORS2 on the VLT. 

The foUowing sections describe the details of the data re- 
duction, calibrations and combination of the different data sets. 

3.1. HST/STIS images 

The HST/STIS images of the host of GRB 990712 were ob- 
tained on April 24, 2000 as a part of a survey of GRB hosts 
(Fruchter et al. 2000b)^. The GRB hosts are observed through 
a clear (unfiltered) aperture (the 50CCD filter, which in the fol- 
lowing analyses will be called the CL filter) and a long pass 
imaging filter, F28 x50LP (called the LP filter). 

The sensitivity of the CL filter extends from 4000 A to 
9000 A, with its peak at --5800 A, which falls within the V 
passband. The sensitivity of the LP filter extends from 5500 A 
to 9000 A, with its peak at ~6000 A. The total integration time 
of the GRB 990712 host was 4080 s in each filters. 

The individual images are combined using DITHER II, a 
package which includes several tasks needed for combining 
dithered HST images. Fruchter & Hook (2002) describe the 
drizzling of WFPC2 images, and we adopt this process for 
the HST/STIS images. This method allows a higher resolu- 
tion in the final images than in the original STIS images. In 
the drizzling of the STIS images of the hosts, the parameters 
pixf rac=0.6 and scale=0.5 were used. This gives an out- 
put pixel size of 0'/0254 pixel" ^ in a 2k x 2k frame. 

3.2. Ground-based optical images 

The DFOSC data in the B, V, R and / filters consist of 1 1, 15, 
15, and 59 frames having total integration times of 12900 s, 
10800 s, 7600 s, and 39950 s, respectively. The individual in- 
tegration times in each filter were not of equal length. 

' http : // www . eso . org/ gen-f ac/pubs/ astclim/ 
forecast /met eo/CIRA/ images /repository/ lossam/ 

^ The FITS files of this and other hosts can be downloaded at 
http : / / www . if a . au . dk/ -^hst / grbjiosts/ intro . html 
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The raw images were bias subtracted after inspection of 
the overscan region of the CCD. Any residual structure in the 
bias level was corrected for using a normalised median filtered 
bias image. The images were flat-fielded using flat fields ob- 
tained from a combination of several twiUght sky observations. 
The selected flat fields produced reduced images in which the 
background signal varied by less than 2%, which is the limit 
for DFOSC data. The reduced images were WCS calibrated 
with a pipeline written by Andreas Jaunsen (ESO, Santiago). 
Approximately 50 reference stars from the USN02 catalogue 
were used to compute the astrometry. From the WCS cali- 
brated images, the shifts and rotations between individual im- 
ages were found. The images were drizzled in much the same 
way as the HST images, not altering the pixel scale so the pixel 
size of the output image is 0'.'39 pixel"^. 

The drizzling method was tested by comparing the final 
drizzled image with a median filtered image. The drizzling 
method gave a higher S/N for faint objects compared to the 
median combined image. 

Another problem to take into account is the fringes that ap- 
pear in the / band images and which depend on background 
level and consequently on the integration time. Although the 
fringes are large scale structures in the DFOSC images and 
should not introduce a large error when combining 59 individ- 
ual images, they can be removed by construction of a fringe 
frame image. Such a frame is produced by replacing the pixel 
values of all stars with a small pixel value, while also replac- 
ing their neighboring pixels and combining the images using 
a threshold rejection. The resulting fringe image was scaled to 
the exposure time of each image and subtracted. This removed 
the fringes from the images in almost all cases. In a few (^5) 
images the background was higher, and the fringes were re- 
moved by subtracting a scaled fringe frame, with a scale factor 
higher than the integration times. 

The photometric reference stars in the field were adopted 
from Sahu et al. (2000) and no re-calibrations of the field in 
the VRI bands were done. Some of the reference stars in Sahu 
et al. (2000) were saturated in the longest DFOSC V, R and 
/ band exposures so the shortest exposures were analysed in 
order to find secondary standard stars in the field which could 
be used as photometric reference stars in the combined frame. 
The magnitudes were derived by performing relative aperture 
photometry on the field with the PHOT package in IRAF. A 
section of a DFOSC image is shown in Fig. 1 where the stars 
denoted A, B, C, and D correspond to the standard stars in Sahu 
et al. (2000), while the numbers represent the secondary stan- 
dard stars. The magnitudes are listed in Table 1. 

For the B band, cahbration images were taken in June 
2001 using tiie Danish 1.5 m with DFOSC. The standard field 
PG 1657 from Landolt (1992) was observed 4 times during the 
photometric night. This field contains 4 standard stars. Only 
B band data were taken of the field containing the host of 
GRB 990712. The uncertainty of the B band calibration is 8% 
as estimated by standard procedures in IRAF. The errors are 
propagated in quadrature assuming that the errors are indepen- 
dent. 

The U band data consist of 4 x 15 min and one 10 min expo- 
sures. The images were reduced in the standard way and coad- 
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Fig. 1. / band field surrounding the GRB 990712 host. North is up and 
east is left as indicated, and the lengths of the arrows are 1'. The stars 
denoted 1-6 are secondary reference stars, and the stars denoted A, B, 
C, and D correspond to the reference stars 1, 2, 3, and 4 in the paper 
of Sahu ct al. (2000). The GRB host itself lies very close to the circle 
surrounding star A, and is not visible in the representation here. The 
magnitudes of the stars are listed in Table 1. 

ded using the drizzle program without changing the pixel sizes. 
The night was photometric according to ESO's monitoring of 
the night sky conditions at La Silla. The zero point, extinction 
coefficient in the U band, and colour term of the night were 
given by Ramana M. Athreya (private communication). The U 
band magnitudes of the reference stars A, B, C, and D were 
measured in the 4 individual images, and the magnitude of the 
host was obtained using relative photometry in the coadded im- 
age using a 2'.'4 radius aperture. This aperture was the same 
used for the DFOSC images. The resulting magnitudes of the 
host in all filters are given in Table 2. The photometric measure- 
ments were checked with SExtractor (Bertin & Amouts 1996), 
and with aperture photometry we found magnitudes consistent 
with those in the Table within la errors. 

3.3. Ground-based near- 1 R images 

For the reduction of the near-IR JHKs SOFT data, a sky im- 
age was constructed from 6-10 object images obtained imme- 
diately before and after each frame. The number of images de- 
pends on the quahty of the sky subtraction which was evaluated 
by eye. The sky subtracted images were divided by a flat field 
obtained from the NTT/SOFI webpages. A flatfield multiplied 
with an illumination corrected image shows variations on the 
order of 0.2% during a month, and the reduction described here 
gave flatfield accuracies of 1%. 

The integration times were 3 1 x 15 s, 40 x 15 s, and 20 x 15 s 
for the H, Ks, and J band respectively. The exposures were di- 
vided into two sets, one for each day, which were reduced sepa- 
rately. Standard star observations were done right after the sci- 
ence exposures. The extinction coefficients were adopted from 
the NTT/SOFI webpages. The corrections for atmospheric ex- 
tinction are small, because the standard star observations were 
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19.84±0.02 


18.76±0.02 


17.48±0.03 
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_ 
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19.16±0.07 




18.21±0.01 


17.87±0.02 


17.48±0.03 
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19.08±0.02 


17.94±0.02 


16.57±0.03 
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18.90±0.02 


18.37±0.02 


17.86±0.03 
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20.72±0.04 


19.73±0.02 


18.68±0.04 
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— 


6 






19.21±0.01 


18.31±0.02 


17.34±0.03 








A 


18.80±0.06 


18.80±0.06 


17.16±0.01 


16.40±0.01 


15.64±0.01 


14.82 ± 0.02 


14.20 ± 0.03 


14.04 ± 0.04 


B 




17.91±0.05 


16.98±0.01 


16.65±0.01 


16.29±0.01 


15.91 ±0.02 


15.60 ± 0.04 


15.57 ± 0.03 


C 




17.69±0.07 


16.46±0.01 


15.97 ±0.01 


15.50±0.01 


14.96 ± 0.02 


14.51 ±0.04 


14.42 ± 0.04 


D 


17.69±0.07 


17.33±0.06 


15.91±0.01 


15.27±0.01 


14.65±0.01 


13.99 ± 0.02 


13.39 ± 0.03 


13.28 ± 0.03 



Table 1. Magnitudes of the secondary standard stars. The "-" signs appear where the magnitudes for the stars have not been derived either 
because they were saturated in some images, which was the case of the B and C stars, or they were faint and therefore gave large photometric 
errors (> 0.08). These were disregarded in the case of the U, B, and near-IR bands of the secondary standards 1-6. The VRI band magnitudes 
of A, B, C, and D are taken from Sahu et al. (2000). 



performed at an airmass different by only 0.1 from the sci- 
ence observations. Likewise, the colour terms in the near-IR are 
small (^0.02). The images of the standard star were analysed 
in order to find the zero point in each filter. The transforma- 
tion equations did not include a colour term, as this produced 
uncertain fits, are given by: Jinst = jl + Jstd, where Jinst is 
the instrumental magnitude and Jgtd is the standard magnitude. 
The transformations for the H and Ks data are similar. The fit- 
ted zero points are : j 1 = 23.06 ± 0.005, hi = 22.937 ± 0.003 
and kl = 22.367 ± 0.004. These zero points agree with the 
values posted on the NTT/SOFI webpages. 

The shifts between the reduced images were found us- 
ing precor, crossdrizz and shiftfind, and the im- 
ages in each field were combined using imcombine. The im- 
ages were combined using a common median zero-point back- 
ground value, and applying a bad-pixel mask obtained from the 
SOFT webpages to reject hot or dead pixels in the combination. 

The magnitudes of the reference stars (the same 4 stars as 
in Sahu et al. (2000)) were measured in 10 individual frames 
for each filter, and the magnitude of the host was found by per- 
forming relative photometry to the 4 stars in the co-added im- 
ages. The magnitudes of the reference stars are presented in 
Table 1, and the near-IR magnitudes of the host are given in 
Table 2 which also presents the fluxes in the various bands. 
The fluxes are obtained by correcting the magnitudes for the 
Galactic reddening Eb-v = 0.033 estimated from the dust 
maps of Schlegel et al. (1998). Then the magnitudes are con- 
verted to AB magnitudes. Due to the faintness of the host in the 
near-IR, the errors are dominated by sky noise. 

The near-IR magnitudes were obtained using the flux en- 
closed within a circular radial aperture of 1'.'2. At optical wave- 
lengths, the flux of the host is not contained within this aperture 
due to worse seeing. In order to get the right colour of the host, 
an offset equal to the difference in the I band magnitude be- 
tween an aperture of 1'.'2 and 2'.'4 was added. This will then 
provide the right I — J colour, as well as a smaUer photometric 
error. The same method was applied for the H and Ks data. As 
long as the colour gradient in the host is negligible, this does 
not change the near-IR magnitudes one would have calculated 
from larger apertures. The colour gradient of the host is negligi- 



fllter magnitude flux (ynJy) 



u 


23.12±0.05 


1.00±0.05 


B 


23.36±0.09 


2.00±0.06 


V 


22.39±0.03 


4.43±0.16 


R 


21.84±0.02 


6.31±0.20 


I 


21.41±0.03 


7.01±0.29 


J 


20.81±0.17 


7.68±1.21 


H 


20.25±0.19 


8.22±1.47 


Ks 


19.98±0.28 


6.74±1.71 



Table 2. Magnitudes and corresponding fluxes of the host of 
GRB 990712 from all ground based observations. The fluxes in 
column 3 are obtained by correcting for a Galactic reddening of 
Eb-v = 0.033, and offsetting to the AB system before converting 
the magnitudes to fluxes. The flux errors do not include the uncer- 
tainty of the Galactic reddening. The R band magnitude is consistent 
with that derived in Hjorth et al. (2000). 



ble at a radius larger than 1" estimated from the morphological 
study of the HST images described in detail in Section 5. We 
find that there are colour gradients in the central 0'.'25, but this 
has no impact on radii larger than 1'.'2. The K band magnitude 
is consistent with the value reported in Le Floc'h et al. (2003). 

4. Spectroscopy 

The host galaxy of GRB 990712 was observed in service mode 
with FORS2 at UT4 of ESO's Very Large Telescope on July 
18, 2001. The exposure time was 10 min, and the grism used 
was GRIS 600zH-23 with order separation filter OG590, giving 
a wavelength range of 7400-10700 A, which includes the red- 
shifted wavelength of the host galaxy's Ha at 9404 A. The slit 
width was set to 1", resulting in a resolving power of approxi- 
mately 1400. The seeing during the observations was around 
(y.'8. The spectrum was reduced in the standard fashion, us- 
ing IRAF. The wavelength calibration was performed using a 
HeNeAr lamp spectrum; the resulting scatter is 0.02 A. 

The resulting spectrum has a bright Ha emission line at 
9404 A which has an observed equivalent width (EW) of 
180±40 A estimated from fitting the continuum level by spline 
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9200 9300 9400 9500 9600 9700 
wavelength (A) 



Fig. 2. Section of the spectrum containing the redshifted Ha line. The 
spectram has been wavelength calibrated, and set to a flux scale cor- 
responding to the continuum level of 7.38 /^Jy estimated from the 
SED analysis. The redshift measured from the Ha line is 0.4331. 
At 9436.4 A a faint emission line from [NH] A6583 is present, cor- 
responding to z = 0.4338 after applying a standard air to vacuum 
correction, and at 9621 A a line from [SII] A6717 corresponding to 
z = 0.4328 can be seen. 



polynomials of different orders. A small part of the spectrum 
is shown in Fig. 2. One clearly sees the Ha line and also at 
low S/N levels, the [Nil] A6583 line and the [SII] A6717. 
The measured EW of the Ha line indicates that the host is a 
young star-forming galaxy. The starburst99 models provide a 
relation between the rest frame Ha EW and the age of a stel- 
lar population (Leitherer et al. 1999). With a rest frame EW of 
125 ±28 A one would expect an instantaneous starburst age of 
^^6 Myr according to the Starburst99 models, assuming solar 
metallicity. A lower metallicity of Z — 0.001 would increase 
this age estimate by a factor of 2. Furthermore, the EW of H/3 
reported in Vreeswijk et al. (2001a) supports a starburst age of 
~6 Myr according to the Starburst99 models. In the case of a 
continuous star formation rate of 1 Mq yr^^ the inferred age 
from the EW is ^60 Myr, i.e. in both scenarios the presence of 
a young population is inferred. 

5. Morphology of the host galaxy 

In the HST images the size of the GRB host is 1'.'3 x0'.'9. The 
STIS images have revealed that the host of GRB 990712 has 
two separate intensity peaks, and neither of these are located 
at the geometrical center of the host. An image and a contour 
plot of the host is shown in Fig. 3. The rightmost knot (south- 
east=SE) is ~1 mag brighter than the left (north-west=NW). 
The GRB occurred in the SE knot within 0'.'048±0'.'080 of the 
center. The afterglow itself will not contribute significantly to 
the flux of the host if a break in the light curve is present around 
one day after the burst as suggested in Bjornsson et al. (2001). 
If no such break occurred, the contamination of the SE knot due 
to the late time afterglow will be ^5%. If a supernova of similar 
brightness as SN 1998bw is present at the time of the observa- 
tions, it would have the magnitude V « 27.0 at z = 0.433. 
This magnitude is calculated assuming a similar late time su- 
pernova light curve contributing to the total flux in addition to 




I ? 3 



Fig. 3. A STIS CL image of the host of GRB 990712 overlayed by 
a contour plot shows that there are two intensity peaks in the image. 
The scale of the plot is 3"x3". The distance between the two peaks 
is 10 drizzled pixels, or 0'.'25, which corresponds to a separation of 
~1500 pc at 2; = 0.433. At the time of observation more than one 
year from the burst the flux contribution from the afterglow will be 
negligible. 



the afterglow. If a SN is present, the contamination of the SE 
knot will be an additional ~5%. 

The colours of the two knots were analysed by perform- 
ing aperture photometry centered on each knot using both the 
CL and LP images. At a redshift of z = 0.433 the difference 
in magnitudes in these filters roughly corresponds to the rest- 
frame B — V colour. The B ~ V colour of the NW knot is 
1.3 ± 0.1, and the colour of the SE knot is 0.9 ±0.1 within an 
aperture of radius 0'.'076 (i.e. 3 drizzled pixels). The colour be- 
comes more red with increasingly larger photometric apertures 
(up to O'.'IS radius, i.e. 6 drizzled pixels), but the colour gra- 
dient in the NW knot is smaller than in the SE knot. At larger 
apertures, the colours will be contaminated by flux from the 
other knot. The differences in colour could be caused by dif- 
ferent ages of two starburst regions, the bright blue knot being 
slightly younger than the fainter. Another explanation could be 
a relatively larger extinction in the NW knot. If the colour ex- 
cess, Eb-v , of the faint knot is larger than 0.34 than for the 
bright knot, this could explain the colour difference. Compared 
to the overall extinction estimated from the SED this is rela- 
tively large. Even though we find no evidence for a large ex- 
tinction from the SED analysis as explained in Section 6, the 
morphology of the host supports the presence of two different 
stellar components. 

Most interestingly, the GRB was coincident with the center 
of the blue SE knot (Bloom et al. 2002b). This location of the 
GRB corresponds to the bluest part of the host, which likely 
links the GRB to a star-forming site. If the two knots were the 
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result of a merging of two components, then one would expect 
to see further evidence of the tidal interaction in, for example, 
luminous tails. This is not apparent in the image, but for low 
mass systems bright tails are probably rare. One should note 
that the surface brightness of tidal tails can be low and thus 
difficult to detect, since the redshift of the host gives a further 
factor (1 + z)^ dimming of the surface brightness assuming a 
standard cosmology. 

6. Spectral energy distribution of tlie liost 

The magnitudes in Table 2 were used to compare to theoretical 
galaxy template spectra from Bruzual & Chariot (1993). This 
was done by the program HyperZ^ described in Bolzonella 
et al. (2000). HyperZ was written mainly for obtaining the pho- 
tometric redshifts of galaxies in large surveys, but it also serves 
the purpose for finding the best matching theoretical galaxy 
template for a given set of broad band observations. The tem- 
plates consist of elliptical, several types of spiral, irregular, and 
starburst spectra at various ages having different star formation 
histories. The time evolution is described by SFR oc exp(-t/T), 
where r is the SFR timescale which increases along the Hubble 
sequence, with r ^ in the case of an instantaneous star- 
burst. The metallicities of the templates are equal to the solar 
value, Z=0.02. Gorosabel et al. (2003a,b) have shown that for 
the GRB 000210 and GRB 000418 hosts, the metalUcity is a 
secondary variable in comparison to the impact of the assumed 
IMF and the extinction law. 

We used both the Miller & Scalo (1979) initial mass func- 
tion (IMF) and the Salpeter (1955) IMF for stellar masses 
between 0.1 and 125 Mq for calculating the templates. The 
Miller & Scalo IMF produces fewer massive stars compared to 
a Salpeter IMF, and at masses below 1 M0 the Miller & Scalo 
IMF is flatter (Miller & Scalo 1979). The largest differences 
between the templates are at the red and near-IR wavelengths. 
Different amounts of internal extinction can also be applied to 
the templates. In this way, the SED fitting allows an estimate 
of the type of galaxy, age, and internal extinction. In the anal- 
ysis the redshift of the templates was fixed to the value of the 
host (z = 0.433). Leaving the photometric redshift Zphot as a 
free variable gives 2;phot = 0.445 ± 0.020, which is consistent 
with the spectroscopic one. This additional free parameter does 
not change any of the resulting values from the best fit besides 
changing the reduced by a small amount. The agreement 
between the photometric and the spectroscopic redshift shows 
that the SED fitting technique is reliable for estimating other 
properties of the host. 

The goodness of the fit is evaluated by the expression: 

2 / ^host,i — k X -Ftomp,i\^ 

where the sum is to be taken over all filters, i. The flux values 
of Fhost are found in Table 2, a"(fhost,i) is the photometric er- 
ror, and A: is a normalization constant. ^temp,i is the flux of the 
template in the filter i, which is calculated using the through- 
put for the given filter and instrument. For the DFOSC and 

^ http : / / webast . ast . obs-mip . f r /hyper z/ 



EFOSC data, the filter transmission was convolved with the 
quantum efficiency of an ESO LORAL CCD in order to calcu- 
late the throughput. For the NTT data, the combined throughput 
of SOFI and the filters were used. 

The SED of the host was best fit by a starburst template. 
This confirms the conclusion in Vreeswijk et al. (2001a), de- 
rived on the basis of emission Unes, that the host is a starburst 
galaxy. The best fit model has a starburst age of 0.255 Gyr 
and an extinction of Ay = 0.15 ± 0.1 using a Salpeter IMF 
for the templates. The error in the extinction was estimated 
from results of the fits for which the per degree of free- 
dom, x^/dof< 2, and all these fits gave an age of 0.255 Gyr. 
Fig. 4 shows the best fits when using templates from a Miller 
& Scalo IMF (thin line) and a Salpeter IMF (thick line) respec- 
tively. Both templates gave the same values for the extinction 
and age for the best fit. The fit to the thin line has a reduced 
xVdof = 2.82 and a fit to the fliick line has x^/dof = 0.959. 
The largest difference between the two templates is in the near- 
IR, where our photometric points have large uncertainties. 

Fitting the SED to other types of galaxy spectra give larger 
values of x^, e.g. x^/dof = 14.3 for the best fit to an irregular 
galaxy template, and y^ldof = 16.2 to an elliptical galaxy. For 
the latter template, all the measured near-IR fluxes were ~3c7 
below the template flux. Generally, all other templates besides 
starburst templates fail to reproduce the flat continuum from 
8000-22000 A, while at the same time fitting the Balmer jump 
at the rest-frame 3646 A. Thus, we infer that the host is most 
likely a starburst galaxy with a stellar distribution similar to a 
Salpeter IMR 

The precision of the age estimation relies upon the accu- 
racy of how well the Balmer jump is sampled. At the redshift 
of the host, this jump will he at ~5200 A. It is seen in Fig. 4 
that with the current set of broad band magnitudes this jump 
is well sampled. Therefore, the age of the dominant population 
of stars is well constrained. However, if more than one popu- 
lation of stars is present in the host, this will change the age 
determination somewhat as explained in section 6.1. 

The extinction found by HyperZ is Ay = 01 5 using the 
extinction curve from Calzetti et al. (2000) appropriate for star- 
burst galaxies. The extinction measured from emission line 
widths was Ay = 3.4^^ 7, which is consistent with a small ex- 
tinction measured from the SED. We also tried to do the SED 
fitting with other extinction curves. The Ay found by HyperZ 
was found not to vary much {Ay <^ 0.2 in all cases) using the 
extinction curve of the Milky Way (MW) from Seaton (1979), 
the Large Magellanic Cloud (LMC) from Fitzpatrick (1986), 
and the SmaU MageUanic Cloud (SMC) from Prevot et al. 
(1984) respectively. Acceptable values of x^/dof< 1.5 for the 
fits were found using the starburst, LMC, and SMC extinction 
curves, so the extinction curve could not be constrained from 
the SED. A larger x^/dof = 2.2 was produced using a MW 
extinction curve, which suggests that the dust in the host is dif- 
ferent than from MW dust. This is in agreement with results ob- 
tained from other GRBs, where the extinction law has been in- 
ferred from studying the afterglows. Jensen et al. (2001), Fynbo 
et al. (2001), Lee et al. (2001), and HoUand et al. (2003) find 
that the SMC extinction law gives a better fit to multiband ob- 
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Fig. 4. Best fit of the GRB host SED to synthetic spectra calculated by the HyperZ program. The synthetic spectra are shown by the two solid 
lines and the dots denote the UBVRIJHKs fluxes from Table 2. The horizontal errorbars indicate the FWHM of each filter The thin line shows 
a template spectrum calculated using a Miller & Scalo IMF while the thick line is calculated using a Salpeter IMF. The values of the photometric 
redshift, age, and extinction are given as inserts in the panel. A starburst extinction curve is used for the calculation of the extinction. 



servations of the afterglows of GRB 000301C, GRB 000926, 
GRB 010222, and GRB 021004, respectively. 



From the SED of the host we can determine the abso- 
lute magnitude in various bands by convolving the restframe 
spectrum of the best fitting template with various filter trans- 
mission curves. This gives AIb = —18.8, My — —19.9, and 
Mfj = —20.3 in the assumed cosmology, corrected for Galactic 
extinction. This absolute magnitude is comparable to that of 
other hosts (Bloom et al. 1998, 2001; Sokolov et al. 2001). The 
luminosity of the galaxy is less than the characteristic lumi- 
nosity L* given by the Schechter function. For field galaxies a 



value of M t 



21 is typically assumed, while for starburst 



galaxies at high redshift this magnitude is highly uncertain, but 
is likely to be brighter than -21 (Lilly et al. 1995). GRB hosts 
have been shown typically to be under-luminous (Le Floc'h 
et al. 2003). From the early optical light curves of the optical 
transient, obtained in the V, R, and / bands when it was still 
bright, Sahu et al. (2000) estimated the luminosity of the under- 
lying host to be of the order of i* . Using high spatial resolution 
HST images Hjorth et al. (2000) refined the luminosity of the 
galaxy to be 0.2L*. 



6. 1. Two population fits 

The scenario described above, where the whole spectrum of 
the galaxy can be represented by one single burst of star for- 
mation is likely too simple. One must expect that more than 
one burst of star formation would contribute to the total mix of 
stars observed. With the knowledge that some GRBs originate 
from collapsing massive stars, we would expect a population 
of younger stars to be present, which is supported by the large 
EW of the Ha line. We therefore investigated whether the SED 
of the host could be explained by a superposition of two popu- 
lations of stars, i.e. a young burst superimposed onto an older 
population. 

The method applied was as follows. We created two new 
SEDs, a red and a blue one, whose sum was the total SED of 
the observed galaxy. The two objects were run though HyperZ, 
finding the best fit templates using the same templates as in the 
single population case. The best fit spectra of the two popula- 
tions were summed and compared to the broad band fluxes of 
the host, and the x^s of the fits were calculated using Eq. (1). 
It was then investigated by iterations if a bluer first population 
plus a redder second population would produce a better fit to 
the observations. 
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This process was done first with two populations of similar 
total flux. In a second run, the observed flux was partitioned 
into 80% for the first population and 20% for the second. 

We found that several two population models were able to 
fit the observations with x^/dof « 1. As shown in the upper 
panel in Fig. 5, the total flux can originate from two rather sim- 
ilar instantaneous burst populations. In this specific case, one 
of the populations has an age of 0.36 Gyr and an extinction of 
Ay = 0.12, while the other has an age of 0.18 Gyr and zero 
extinction. The fit of the summed spectrum to the observations 
is x^/dof=0.66. 

Another scenario could also explain the properties of the 
host. The lower panel in Fig. 5 shows the result of the first pop- 
ulation of stars being a 0.36 Gyr old starburst, while the flatter 
spectrum corresponds to a less luminous second population, 
which is a 52 Myr starburst. Both these spectral templates have 
Ay = 0.00. The fit of the summed spectra to the observations 
gives x^/dof = 0.75. We therefore conclude that if a second 
population of stars has a significant contribution to the total 
flux, the age will be 50-200 Myr, found from acceptable val- 
ues of the x^/dof fits. Considering the large observed Ha EW 
we find that a young population (~50 Myr) is preferable. Even 
younger ages for the second population can not be ruled out. 
By constructing templates of very young populations of stars, 
we can estimate the total flux allowed from such a population. 
The constraining factor is the weak blue continuum observed. 
If a stellar population with an age of lOMyr is present, the total 
(bolometric) flux it emits is less than 5% of the total observed 
flux from the host. 

7. Star-formation rate 

The continuum in the UV part of the spectrum (1500-2800 A) 
mainly comes from young OB stars, and a relation between the 
UV flux and the SFR can be derived by comparison of observed 
spectra to synthetic model spectra (K98). 

We interpolate between the observed flux in the U and the 
B band, and the flux at the wavelength 2800(1 + z) A is esti- 
mated assuming a powerlaw spectrum, Fhost oc i''^ . We do not 
consider any best template fit for the calculation of the flux at 
this wavelength. The total luminosity at 2800(1 + z) A can be 
calculated given the cosmological model, and converted to an 
overall SFR of the host. K98 gives the relation between the UV 
luminosity and the SFR for a Salpeter IMF. 

SFR(Mq yr-1) = 1.4 x lO^^s ^ l^^^ [erg s^^ Hz-i](2) 

This expression is adequate for galaxies with constant stel- 
lar formation over time scales of at least 10® years. The coef- 
ficient that finks the SFR with the UV luminosity (1.4x 10"^* 
in the case of Eq. (2)) is significantly lower in younger stel- 
lar populations present in starburst regions. Thus, this SFR has 
to be considered as an upper limit. However, in starburst re- 
gions dust is expected to be present, and a correction for dust 
extinction is necessary in order to find the UV flux. We have 
in the SED analysis found only a small global extinction ef- 
fect, but it could be locally higher if the medium is clumpy. 
For instance, the appreciable extinction inferred along the line 




5.0'10^ 1.0*10'^ 1.5'10'* 2.0*1 0"^ 
observed wavelength (A) 







^ if i 1 

- Jl 







S.O'IO^ 1.0*1 0'* I.S'lO'* 2.0*10'^ 

observed wavelength (A) 



Fig. 5. Two examples of the best fits for two populations. In both pan- 
els the best fit spectra for each of the components are shown along with 
the summed spectrum. In the upper panel, the lower spectrum corre- 
sponds to a 0. 1 8 Gyr starburst population, the middle spectrum to a 
0.36 Gyr starburst population, and the upper spectrum to a sum of the 
two. The fit has x^/dof = 0.66. In the lower panel, the lower spectrum 
corresponds to a 52 Myr starburst population, the middle spectrum to 
a 0.36 Gyr starburst population, and again the upper spectrum to a sum 
of the two. This fit has x^/dof = 0.75. 

of sight to some afterglows derived from the curvature of the 
near-IR/optical afterglow SEDs in e.g. GRB 000301C in Jensen 
et al. (2001), and GRB 000926 in Fynbo et al. (2001) seem to 
indicate that GRB progenitors tend to be embedded in dusty 
local regions. High Nh values derived from X-ray afterglow 
spectra also supports this scenario (Piro et al. 2002). This will 
lead to an underestimate of the local extinction, and hence of 
the SFR. This opposite effect may mitigate the SFR overesti- 
mation entailed in Eq. (2). The estimated SFR value has to be 
considered as the sum of the SFR contributions coming from 
the two potential stellar populations. 

We calculate a flux at the rest frame 2800 A 
of 1.65±0.04 fiJy which translates into a SFR of 
1.3±0.3 M© yr"\ which is similar to the SFR of a 
galaxy such as the Milky Way. The errors due to the flux 
measurements and interpolation between the two bands are 
insignificant compared to the uncertainty of the conversion 
factor, which is ^30%. CoiTecting for the extinction of 
Ay ~ 0.15 using the starburst extinction curve from Calzetti 
et al. (2000) gives a slightly larger flux at 2800 A in the rest 
frame, and the SFR is a bit larger: 1.6±0.3 M© yr~^. 
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We also estimate the SFR from the Ha line flux. K98 gives 
the relation between the Ha luminosity and the SFR: 

SFR (Mq yr-i) = 7.9 x lO"*' x [erg s'^] (3) 

We can find the line flux, /iine> and the luminosity, ^Kne, once 
the continuum flux /cont at the wavelength of the redshifted Ha 
line A, and the equivalent width EW is known. 

/line = /cont X c/A^ X 10"^° X EW 

hine = /line47rd|, (4) 

Since the spectrum has not been flux calibrated, we esti- 
mate the continuum level by interpolation between the / and 
J band fluxes assuming a powerlaw spectrum. We find a con- 
tinuum level of (7.38 ± 0.44) /jJy. Using Eq. (4) we find 
/line = (4.50 ± 1.00) X lO"^*" erg s~^ cm"^, which gives a 
SFR=(2.8±0.7) Mq yr-^ The error includes the uncertainty 
of the conversion from flux to SFR, which is 30%. The two 
calculated SFRs are consistent with each other within 2a. 

The larger SFR inferred from the Ha line compared 
to the SFR found from the UV flux could be due to 
dust extinction, which is stronger in the UV region. The 
ratio between the Ha line flux and the H/3 line flux, 
(1.33±0.20) X 10^1^ erg cm^^ ^-i jjj Vreeswijket al. (2001a), 
is 3.38±0.90, while the expected ratio in HH regions in the 
case of no extinction is 2.85 (Osterbrock 1989). The ratio be- 
tween the observed line fluxes thus corresponds to a magni- 
tude difference of 0.19 ± 0.31. Using the extinction curve from 
Calzetti et al. (2000) an extinction of Ay = 0.60 ± 0.99 is in- 
ferred from the line ratio. The calculated extinction is the same 
when using the MW extinction curve from Fitzpatrick (1999), 
while in the case of an SMC extinction curve one would find 
Av = 0.41 ± 0.60. The extinction is therefore consistent with 
the small value indicated by the SED analysis. 

The SFR of this host is in the same range as the SFRs 
found for other hosts through their rest frame UV flux, which 
typically gives SFRs < 10 Mq yr"^ (Fruchter et al. 1999; 
Bloom et al. 1998; Djorgovski et al. 2001). The largest SFRs 
found from optical methods to date are 20 Mq yr^^ for the 
GRB 990703 host (Djorgovski et al. 1998) and 55 M© yr"! for 
the GRB 000418 host (Bloom et al. 2002a), although a smaller 
SFR based on the UV region of the latter host has been inferred 
(Gorosabel et al. 2003b). It must be pointed out that these SFRs 
are strictly lower limits to the true SFRs since the reported val- 
ues are not corrected for extinction by dust in the hosts. Radio 
and sub-mm data suggest SFRs one or two orders of magnitude 
larger than the optical inferred SFRs for a sample of GRB hosts 
(Berger et al. 2003). 

8. Discussion and conclusions 

From broad band magnitudes in UBVRIJHKs filters we have 
examined the SED of flie host of GRB 990712. Comparing 
this SED with model templates of different galaxy types, we 
found that the host is a starburst galaxy with an extinction of 
Av = 0.15. With spectroscopic observation of the host we 



calculated the extinction Ay = 0.6 ± 0.99 from the Ha/H/3 
Une ratio, confirming a small extinction value. 

In the collapsar scenario the progenitor of the GRB may be 
embedded in a molecular cloud having a much larger extinction 
due to the surrounding dust. Thus, even though we can estimate 
the overall internal extinction in the host in the particular case 
of GRB 990712, it is not possible to say anything about the ex- 
tinction in the line of sight towards the burst itself. It could well 
be much higher. However, analyses of several afterglows have 
failed to reveal a very high extinction {Ay > 1). According 
to Galama & Wijers (2001), the expected visual extinction of 
the afterglows should be much higher when compared to the 
column density inferred from the X-ray afterglows of several 
bursts. They argued that dust can be destroyed along the Une 
of sight towards the burst making the visual extinction appear 
smaller. Waxman & Draine (2000) have calculated that dust 
can be destroyed out to a distance of 10 pc from the burst site. 
The extinction inferred from the SED fitting is an overall aver- 
age extinction of the entire galaxy. Therefore it is necessary to 
investigate further the relation between the smaU extinction in- 
ferred from the optical light curves of some GRBs, Ay < 0.2 
(Andersen et al. 2000; Galama & Wijers 2001; Jensen et al. 
2001; Fynbo et al. 2001; Stanek et al. 2001), and the extinction 
in the host itself. 

The SED of the host is similar to that of a starburst popu- 
lation with an age of 0.26 Gyr at a redshift of 0.43. This age 
is still consistent with a merging neutron star scenario as the 
progenitor of the GRB. It is now known that some of the long- 
duration GRBs are associated with coUapsing massive stars 
(Stanek et al. 2003; Hjorth et al. 2003). Considering that the 
life times of the most massive stars are of the order of a few 
Myr, a small age of the star burst is expected. It was therefore 
investigated whether two distinct populations were able to fit 
the broad band observations of the host. It was found that in 
such case, the best fits were produced by a younger starburst 
population with an age of 50-180 Myr, and zero extinction. 
We consider the lower limit to be more likely given the large 
Ha EW found in the spectrum, which suggests presence of a 
young stellar population with an age of 6-60 Myr depending 
on the star formation history. 

From the analysis of the SED we found that a Salpeter IMF 
was able to reproduce galaxy spectral templates corresponding 
to the observed fluxes. 

We calculate the SFR by estimating the rest frame flux at 
2800 A. The SFR=1.3±0.30 M© yr^^ is not very large, and 
correcting for internal extinction in the host does not increase 
the SFR much. This relatively small SFR is comparable to that 
of other GRB hosts found from using the same UV-SFR esti- 
mator. Considering that the host is less luminous than an M* 
galaxy, this SFR is relatively high compared to present day 
galaxies. As the host is a 0.2L* galaxy we find a SFR per L/L* 
of~5 Mq yr^i {L/L*)-^. 

Comparison wifli flie SFR of 2.8±0.7 Mq yr"! found 
from the Ha line flux implies that there may be moderate ex- 
tinction present in the host. The dust may be distributed in a 
clumpy medium, where most of the UV flux is absorbed. 

The analysis of the morphology of the host showed that it 
has two knots of different colours. This colour difference could 
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be due to two bursts of star formation. This interpretation is 
consistent with the large Bahner break in the SED which sug- 
gests the presence of an older population together with the large 
Ha EW suggesting the presence of a very young population of 
stars. 

Most importantly, the location of the burst was in the bluest 
part of the host galaxy, which supports the recent observations 
that the long-duration GRBs are Unked to sites of formation of 
massive stars. 
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